Momentum Conserving Model with Anomalous Thermal Conductivity in Low 

Dimensional Systems 
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Anomalous large thermal conductivity has been observed numerically and experimentally in one- 
and two-dimensional systems. There is an open debate about the role of conservation of momentum. 
We introduce a model whose thermal conductivity diverges in dimension 1 and 2 if momentum is 
conserved, while it remains finite in dimension d > 3. We consider a system of harmonic oscillators 
perturbed by a non-linear stochastic dynamics conserving momentum and energy. We compute 
explicitly the time correlation function of the energy current Cj(t), and we find that it behaves, for 
large time, like t~ d ^ 2 in the unpinned cases, and like ^ _<i / 2_1 when an on site harmonic potential 
is present. This result clarifies the role of conservation of momentum in the anomalous thermal 
conductivity in low dimensions. 
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When a small gradient of temperature VT is applied 
to a material, we expect that, in the steady state, the 
heat current satisfies Fourier's law 



< J > = 



-kVT 



where k is the conductivity of the material. There has 
been interest in the question of validity of Fourier's law 
for low dimensional systems (d < 2), where standard ar- 
guments based on linear response break down (cf. 0, Q 
for a review on the subject). Anomalous large conductiv- 
ity is observed experimentally in carbon nanotubes and 
numerically in Fermi-Pasta-Ulam (FPU) systems with- 
out pinning (on-site potential), where numerical evidence 
shows a conductivity diverging like N a , with a < 1 in 
one dimension, and like log N in dimension 2 (cf. Q and 
references therein) . If some non-linearity is present in the 
interaction, finite conductivity is observed numerically in 
all pinned cases, and it is expected in all 3 dimensional 
cases P, • Consequently it has been suggested that con- 
servation of momentum is an important ingredient for the 
anomalous conductivity in low dimension 3J. 

In insulating crystals heat is transported by lattice vi- 
brations, and since the pioneering work of Debye, systems 
of coupled anharmonic oscillators have been used as mi- 
croscopic models for heat conduction. Non-linearity is 
extremely important. In fact, in the linear case the aver- 
age energy current < J > is independent of the length N 
of the system, i.e. the conductivity kn diverges like A?" 0. 
In fact in the harmonic crystal the normal modes of the 
vibrations (phonons) do not interact and follow ballistic 
motion. A finite asymptotic conductivity instead should 
result eventually by the diffusive behavior of phonons due 
to phonon-phonon interaction caused by anharmonicity. 

Since conductivity in non-linear systems is difficult to 
compute or estimate analytically, it is natural to model 



the nonlinearities by stochastic perturbations of the lin- 
ear dynamics. In some sense these stochastic perturba- 
tions simulate (qualitatively) the long time (chaotic) ef- 
fect of the deterministic non-linear model. 

We study in this Letter a stochastic model where con- 
ductivity can be explicitly computed, and diverges in di- 
mension 1 and 2 when momentum is conserved, while 
it remains finite in dimension 3. So this is the only ex- 
plicitely solvable model that has a behavior qualitatively 
consistent with numerical simulations. 

We consider a system of harmonic (linear) coupled os- 
cillators where the Hamiltonian dynamics are perturbed 
by a random exchange of momentum between nearest 
neighbor atoms. The random exchange of momentum 
conserves total momentum and total energy. We con- 
struct this noise with a diffusion on the surface of con- 
stant kinetic energy and momentum. Because of the con- 
servation laws, this noise introduces a certain nonlinear- 
ity in the model. 

We compute explicitly the time correlation function at 
equilibrium of the energy currents Cj(t) (cf. equation 

and we find that, as t hoo, Cj{t) ~t~ d / 2 if 

the system is unpinned, while Cj(t) ~i _d / 2_1 if an on- 
site potential is present. Conductivity, defined by Green- 
Kubo formula, is then finite only in dimension d > 3 
or for the pinned system. This indicates a divergence 
of the conductivity of the finite system kn as N 1 / 2 in 
the unpinned one-dimensional case, and as log N in the 
unpinned two-dimensional case. 

Other explicitly solvable models have been proposed 
before as perturbation of the harmonic chain (in 0, Q 
only the number of particles is conserved, and in Q only 
energy and the number of particles are conserved). In all 
these models, conductivity is always finite. 

In order to compute the conductivity by the Green- 



2 



Kubo formula, we consider the dynamics of the closed 
system of length N with periodic boundary conditions. 
The Hamiltonian is given by 

Hn = 2 S t p - + qx ' ^ 7 - aA )q-] • 

X 

The atoms are labeled by x £ Tfy, the d-dimensional 
discrete torus of length N. We denote with V and A 
respectively the discrete gradient and the discrete Lapla- 
cian on Tfj. {q x } are the displacements of the atoms 
from their equilibrium positions. The parameter a > 
is the strength of the interparticles springs, and v > is 
the strength of the pinning (on-site potential). 

We consider stochastic dynamics where the probability 
density distribution on the phase space at time t, denoted 
by P(t,qvp), evolves following the Fokker-Planck equa- 
tion (cf. U) 



dP 



VP 



Here L = A + jS is the generator of the process, L* the 
adjoint operator, A is the usual Hamiltonian vector field 

A = {Px ■ 9 qx - \{vl - aA)q x ] • d Px } 



while S is the generator of the stochastic perturbation 
and 7 > is a positive parameter that regulates its 
strength. The operator S acts only on the momentums 
{p x } and generates a diffusion on the surface of constant 
kinetic energy and constant momentum. This is defined 
as follows. For every nearest neighbor atoms x and z, 
consider the (d — l)-dimensional surface of constant ki- 
netic energy and momentum 

§e, P = {(Px, Pz) £ R 2d : p x + p z = 2e; p x + p z = p} . 
The following vector fields are tangent to § e ,p 

xti% = (pL-K)(d Pi - d Pi ) - ( P i-pi)(d pi - d pi ). 

so Yli j— i (-^x,z) 2 generates a diffusion on S e , p . In d > 2 
we define 



1 d 

s =2(d~ij^^( x ^ 



i,j,k 



where ei, . . . ,e<j is the canonical basis of Z d . Observe 
that this noise conserves the total momentum J2 X Px and 
energy Hn, i.e. 



S ^ p x = , SHn=0 



In dimension 1, in order to conserve total momentum 
and total kinetic energy, we have to consider a random 



exchange of momentum between three consecutive atoms, 
and we define S = | X^eT 1 O^x) 2 , where 

Y x = (p x ~Vx+i)dp !C _ 1 + (p x+ i -Px-i)d Px 
+ (Px-i ~Px)d Ptc+1 

which is vector field tangent to the surface of constant 
energy and momentum of the three particles involved. 

These dynamics can also be written in terms of the 
solutions of the stochastic differential equations 

dp x = - (W - A) q x dt + 27Ap x dt + y/j dn x (t) (1) 

where of course q x = p x and n x (t) are defined by the 
Ito's stochastic integrals 

nx(t)= 27f=T £ E / (^>x) ( S ) d< y (.) 

V ||y-x||=l i,j 

Here w^ y (t) = u>^ x (i) are independent standard Wiener 
processes. In d = 1 the expression is similar with the 
term 27Ap x replaced by ('j/6)A(4p x +p x +i +Px-i)- 
Defining the energy of the atom x as 

1 2 a ^-^ 9 

e x = ^Px + 4 2^ ( q y ~ qx )" + 2 qx ' 

y:|y-x| = l 

the energy conservation law can be read locally as 



ex(t) - e x (0) - £ ( J — e.,x(t) - Jx,xW*)) 
fc=i 

where J x , x +e fc (t) is the total energy current between x 
and x + efc up to time t. This can be written as 

Jx.x+e, (*) = / jx,x+e* (s) ds + Af x ,x+e fc (*) • (2) 
JO 

In the above M X}X+ek (t) is the Ito's stochastic integral 
defined by 

M x , x+efc (i) = ^L^jl (^ +et e x ) (*) dw%{s) 

The instantaneous energy currents j x , x + efc satisfy the 
equation 

d 

Le x = ^ ] (jx-e t ,x — Jx,x+e fc ) 
fc=l 

and can be written as 



Jx,x+e fe — Jx.x+ei. + 7.7 



fc 1 (Jx,x+e fc 



(3) 



The first term in © is the Hamiltonian contribution to 
the energy current 



OL 

Jx'.x+e, = - "jjtox+efc - fix) ' (Px+e fc + Px) (4) 
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while the noise contribution in d > 2 is 



13 



x,x+e fc 



-7Ve fc P x 



(5) 



and in d = 1 is 



V(Px-l,Px,Px+l) =^[Px+l +4pl+pl_ 1 +p x+1 p x _i 

- 2p x+1 p x - 2p x p x _i] 

Consider the dynamics of the closed system on in 
microcanonical equilibrium. The microcanonical distri- 
bution is usually defined as the uniform measure on the 
energy surface Ti — N d e, for a given e > 0. Our dynam- 
ics conserve also (EPx) 2 + KE^lx) 2 - Notice that the 
dynamics is invariant under the change of coordinates 
Px = Px E y Py and q x = qx E y %■ Consequently, 
without any loss of generality, we can fix E Px = and 
E Qx = in the microcanonical measure. 



Let us define 3 ei = Exi*,x+ei = E x ix,x+ei- We are 
interested in computing the correlation function: 



C 1A (t) = \in^—E(3 ei m ei (0)) (6) 

where E is the expectation starting with the microcanoni- 
cal distribution defined above. By explicit calculation we 
can solve the equation 

(A - Ly^ ei = --^ ffA:A r(x- y)p x • q y (7) 

x,y 

where g\^(x) is the solution of the equation 



2A 

— Sa .jv(x) - 4Ag A ,jv(x) = (<J(x + ei) - S(x - ei)) (8) 

7 



on for d > 2, while in one dimension it solves 



J 



2A 1 

— 9x.n(x) - -A [4 5AjA r(x) + g x . N (x + 1) + g x . N (x - 1)] = (S(x + 1) - 5(i - 1)). 

7 " 



r 



(9) 



Consequently, for A > 0, we can write the Laplace trans- 
form of Ci t i(t) as 



e- Xt C 1A (t) dt = lim U ei (A - i) _1 a ei ) w (10) 



where < • • • >jv denotes the expectation with respect to 
the microcanonical measure. 

Substituting Q in (|1Q [I and using equivalence of en- 
sembles estimates (see |2| for the details) we have 



e" At C M (i) dt 



2dj 



^3a(z) [r(0,z + e 1 )-r(0, Z -e 1 )] 



(11) 



where L is the kernel of the operator (vl — qA)' 1 on 
Z d , while g\ is the solution of equations {HJ in Z d or 
in Z. We compute explicitly (| 1 1 |t and invert the Laplace 
trasform obtaining 

CnW = t4j / (5 fc ^(k)) 2 dk (12) 

47r « y[o,i]" 

where w(k) = (f + 4aEj=i sin 2 ^^ )) 1 / 2 is the disper- 
sion relation of the system, and 



^ = f 8E"=iSm 2 (7T^) ifd>2 
^ V ; \ 4/3sin 2 (7rfc)(l + 2cos 2 (7rfc)) if d = 1 



' (13) 



Since around k = 0we have ^(k) ~ k 2 and 

(9feiw(k)) 2 ~ (a/c 1 ) 2 (^ + a 2 k 2 ) -1 , we have the following 

asymptotic behavior 



C w (t) 



t -d/2 if j/ = , 



t^oo | t-d/2-1 if v > 



(14) 



By Green-Kubo formula the conductivity in the 
direction ei is given by 

k 1 ' 1 = , lim J im E ( Jx - x + e i W J °- e i W) (15) 

t^oo 2e I N— >cx) * — ' v ' 



By explicit calculation one can show (see |8j for the 
details) 



,1 7 . 1 



, , , , C ltl (t) dt 
d 2e z J 



7 , 1 



d 87r 2 c?7 7[o,i] d V'C 4 ) 



(16) 



dk. 



By (|14p. if <i > 3 or if ^ > 0, the integral on the right 
hand side of (|16l) is finite. 

If v = and d < 2, by 1|14|) the time integral 
in l|16fl diverge and conductivity is infinite. Follow- 
ing Ref. (page 46), one can estimate the depen- 
dence of the conductivity of the finite open system of 
length N with thermic baths at the boundary imposing 
a temperature gradient, by restricting the time integral 
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in l|16|) to times smaller than the "transit time" N/v s , 
where v s is the sound velocity of the lattice defined as 
v s = liming d k iu(k) — 1. This gives a finite size con- 
ductivity kn diverging like TV 1 / 2 in dimension 1, and like 
logiV in dimension 2. 

Discussion. The exact results presented in this letter 
concerning the stochastic model we introduced give some 
indications about the role of conservation of momentum 
and of confinement (pinning) in heat conduction for the 
nonlinear deterministic Hamiltonian case. In fact the de- 
cay of the energy current correlations and consequently 
the behavior of the conductivity that we proved in our 
stochastic model are (qualitatively) the same as those 
indicated by numerical simulation for the deterministic 
non-linear FPU models. Furthermore a recent paper [l(Tj | 
on the one dimensional unpinned purely quartic FPU 
model suggest the same decay of the time-correlations 
of energy current as in our stochastic model. 

In the one-dimensional case we can give the following 
euristic explanation of the effect of the noise in these 
harmonic systems. In deterministic harmonic systems 
the energy of each mode is conserved, in both pinned 
or unpinned chain; so if modes are created by initial or 
boundary conditions, they cannot interact and they move 
ballistically. This causes ballistic transport and diverg- 
ing conductivity, in both cases (cf. IlLUl). The effect 
of the energy-momentum conservative noise we have in- 
troduced is to scatter modes randomly with rate propor- 
tional to k 2 , for small wave number k. The velocity of the 
k mode is given by the gradient of the dispersion function 
Vw(fc). In the unpinned chain Vw(fc) ~ 1 for small k, so 
small wave number modes have little probability to be 
scattered, and their movement results in a ballistic con- 
tribution to energy transport, while modes with large k 
scatter fast and consequently they diffuse. Properly av- 
eraging over all modes one obtains a current proportional 
to TV -1 / 2 , i.e. a conductivity diverging like jV 1 / 2 . In the 
pinned chain small wavenumber modes move very slowly 
[Vu;(fc) ~ k], so there is a high probability they will be 
scattered and then diffuse while they cross the system. 
Consequently in this case conductivity is finite. 

In |2j we considered the unpinned 1-dimcnsional har- 
monic chain with noise that conserves only energy, and 
prove that conductivity is finite in any dimension. In 
this last case all modes are scattered with constant rate 
(independent of k). 

In non-linear FPU type of interaction, a behavior 
Kpf ~ N a , with < a < 1 is observed numerically. 



But numerical simulations are not conclusive about the 
value of a and there is an intense debate in the litera- 
ture on this value (cf. @). As suggested recently by Livi 
|l2l | , in one-dimensional systems the value of a may de- 
pend on the specific non-linearity of the interaction (un- 
like the logarithmic behavior of the two-dimensional sys- 
tems). The non-linearity creates some scattering of the 
longwave modes, which results in a breaking of the bal- 
listic transport, and in a superdiffusive behavior of these 
modes. An extreme case is given by the 1-dimensional 
coupled-rotors model, which is an example of a non-linear 
chain that conserves momentum and has finite conduc- 
tivity :13||. In this example, rotobreathers (isolated ro- 
tors with high kinetic energy that turn very fast) scatter 
waves that try to pass through them (cf. 
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